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Abstract 

The complex Hilbert space of standard quantum mechanics may be 
treated as a real Hilbert space. The pure states of the complex theory 
become mixed states in the real formulation. It is then possible to gen- 
eralize standard quantum mechanics, keeping the same set of physical 
states, but admitting more general observables. The standard time re- 
versal operator involves complex conjugation, in this sense it goes beyond 
the complex theory and may serve as an example to motivate the general- 
ization. Another example is unconventional canonical quantization such 
that the harmonic oscillator of angular frequency uj has any given finite 
or infinite set of discrete energy eigenvalues, limited below by hco/2. 

1 Introduction 

There are well known mathematical arguments saying that the Hilbert space of 
quantum mechanics could be real or quaternionic, as alternatives to the standard 
complex theory . The real case was dealt with mainly by Stueckelberg and 
collaborators, who concluded that it is essentially equivalent to the complex 
|To}l (see also JH], 0). The interest in the quaternionic case is more alive 
A more exotic subject is octonionic quantum theory [|l7| , [isf . 
The compromise proposed here is a genuine extension of the complex theory, 
but is not quite the full quantum theory on a real Hilbert space. The set of 
physical states is taken to be exactly the same as in the complex theory, but the 
complex Hilbert space is reinterpreted as a real space, and the set of observables 
is enlarged from the set of all complex Hermitean matrices to the set of all real 
symmetric matrices. 

There exists some physical motivation for such a generalization in the fact 
that the time reversal operator T is antilinear in the complex theory. All trans- 
formations involving time reversal are antilinear, among them the fundamental 
CPT symmetry of quantum field theory. It is true that the effect of time re- 
versal can be described easily enough in standard quantum theory, but strictly 
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speaking, as soon as time reversal is introduced, the step from the complex to 
the real Hilbert space has already been taken j| . 

As another example, it is shown in Section |^ below how to represent the 
canonical Poisson bracket relation {x,p} = 1 in terms of operators on a finite 
dimensional real Hilbert space. As is well known, the canonical commutation 
relation [x,p] = —ihl can not be represented on a finite dimensional complex 
Hilbert space, simply because the commutator on the left hand side must then 
have zero trace, while the identity operator / on the right hand side has nonzero 
trace. The argument does not apply in the real Hilbert space, because there the 
operator J = il has zero trace. 

The main argument of Stueckelberg for the equivalence of real and complex 
quantum mechanics is the need for an uncertainty principle. On the real Hilbert 
space it is very useful, if not strictly necessary, to have an operator J commut- 
ing with all observables and having the property that J 2 = — 1, if one wants to 
derive a general inequality for the product of the variances of any pair of observ- 
ables. However, the argument is not compelling, partly because, as Stueckelberg 
points out, there might in principle be one separate operator J for every pair of 
observables, and partly because quantum mechanics makes sense even without 
a general uncertainty principle. In the example of Section || below, the uncer- 
tainty principle for position and momentum holds indeed in all physical states, 
even though the position and momentum operators both anticommute with the 
operator J defining the complex structure. 

Stueckelberg and collaborators also discussed field quantization with fields 
that are either linear or antilinear, in the sense of commuting or anticommuting 
with J. Again their main conclusion is that quantization with antilinear fields is 
impossible for bosons, and possible for fermions but then essentially equivalent 
to quantization with linear fields, so that the real case reduces to the complex 
case. If the conclusion is valid in the present case, it means that the unconven- 
tional quantization of the harmonic oscillator does not lead to any interesting 
new quantum field theory. However, one should perhaps reexamine the argu- 
ments, keeping in mind in particular that there might be several different square 
roots of —1, as discussed briefly in Section ||. 

If antilinear field operators are ever going to be useful, it would most likely 
be in the quantization of the Dirac field. If the Dirac matrices are chosen real, 
then the Dirac equation is seen to be a real equation, and it is not unnatural to 
go one step further and formulate the quantum field theory in terms of a real 
Hilbert space, with linear or (possibly?) antilinear field operators, symmetric 
with respect to the real scalar product. The standard complex notation, both 
for the field equation and for the Hilbert space, hides the fact that the theory 
contains several i = \J — 1 that are logically different, although they are identical 
by the notation. 

One i is the generator of electromagnetic gauge transformations in the Dirac 
equation for a charged particle, it commutes with the mass term in the equa- 
tion. A second i appears in the massless Weyl equation, it generates chiral 
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gauge transformations, and does not commute with the Majorana mass term. 
For this reason, the Majorana and Dirac mass terms are claimed to be differ- 
ent, although one might just as naturally have concluded that there are two 
different i's and only one kind of mass term. A third i turns up in the Fourier 
transformation connecting the position and momentum representations of the 
fields, by definition it commutes with the field operators. The fourth i, acting 
on the complex Hilbcrt space where all the fields act as operators, need not in 
principle be identified with any of the three. 



2 Complex and real Hilbert spaces 

For simplicity, we will consider here mostly finite dimensional Hilbert spaces. 
The complex Hilbert space C D of complex dimension D corresponds to the real 
Hilbert space R 2D of real dimension 2D. The imaginary unit i on C D is then 
a linear operator J on JH 2D with J 2 = —I. Here / is the identity operator on 
R 2D . We define the correspondence so that we have for D = 2, as an example, 
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Then J is an antisymmetric matrix, 
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has a real part which is the symmetric scalar product T ?/> = ip T (l> on R 2D , 
and an imaginary part which is the antisymmetric symplectic scalar product 
-(j, T Jrp = ^ T J(j) on R 2D . 

To any complex D x D matrix A corresponds a real 2D x 2D matrix, which 
we choose to call by the same name A. The correspondence is such that, for 
example, 
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In tensor product notation we may write 



?Mt tn; »Ht z)- <*> 

The Hermitean conjugate of the complex matrix A corresponds to the trans- 
posed A T of the real matrix A. The distinguishing property of those real ma- 
trices that correspond to complex matrices, is that they commute with J. Any 
real 2D x 2D matrix A can be written in a unique way as A = A + + A- , where 
A + J = JA + and A_ J = —JA_, in fact the explicit solution is 

A± = \ (A T JAJ) . (6) 

A + is complex linear. A~ is complex antilinear, thus it is a product of complex 
conjugation and a complex linear operator. In the AD 2 dimensional space of all 
real matrices, the complex linear and the complex antilinear matrices form two 
complementary subspaces of complex dimension D 2 and real dimension 2D 2 . 

The complex D x D matrix A is Hermitean if A^ = A and unitary if A^ = 
A^ 1 . The real 2D x 2D matrix A is symmetric if A T = A, antisymmetric if A T = 
—A, orthogonal if A T = A^ 1 , and symplectic if A T J = JA^ 1 . An orthogonal 
matrix is symplectic if and only if it commutes with J. Thus, the complex 
Hermitean matrices correspond to those real matrices that are symmetric and 
commute with </, whereas the complex unitary matrices correspond to precisely 
those real matrices that are orthogonal and symplectic. 

In other words, an orthogonal matrix is a real linear operator on H 2D which 
is invertible and preserves the ordinary real scalar product <j> T ip. A symplectic 
matrix is invertible and preserves the symplectic scalar product —<f> T Jip. And 
a unitary matrix is a complex linear operator on C D which is invertible and 
preserves the complex scalar product cf^ip. (In the finite dimensional case, but 
not in the infinite dimensional case, the invertibility is a consequence of the 
preservation of scalar products.) 

An infinitesimal linear transformation U = I + eA on R 2D , with e infinitesi- 
mal, is orthogonal if and only if the generator A is antisymmetric, A T = —A. It 
is symplectic if and only if A T J = — J A, which means that the matrix B = J A 
is symmetric. Equivalcntly, A = — JB with B symmetric. U is both orthogonal 
and symplectic if and only if A = — JB = —BJ with B symmetric. 

The dimension of the orthogonal group 0(2D) is 2D 2 — D, the dimension 
of the symplectic group Sp(2D) is 2D 2 + D, and the dimension of the unitary 
group U(D) is D 2 . 
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3 Quantum Mechanics 



3.1 Observables and probabilities 

In standard quantum mechanics a pure state of a given physical system is rep- 
resented by a unit vector ip G C D , or equivalently by the Hermitean density 
matrix p = ipip^, which is a projection operator, since p 2 = p. An observable 
of the system is represented by a Hermitean D x D matrix A, and the theory 
predicts the expectation value in the pure state ip as 

(A) = ^Atp = Tr(pA) . (7) 

More generally, any pure or mixed state is represented by a Hermitean density 
matrix p which is positive definite, i.e. has nonnegative eigenvalues, and has 
unit trace, Trp = 1. The expectation value of the observable A in this state is 
(A) = Tr(pA). The theory also predicts the variance of A, var(A) = (AA) 2 , as 

^(A) = ({A-(A)f) = {A 2 )-{A) 2 . (8) 

A has a sharp value in the state p if and only if var(A) = 0. 

This probability interpretation makes sense because of the spectral theorem 
for Hermitean matrices, which guarantees the existence of the spectral repre- 
sentation 

N 

A = J2 a nPn • (9) 

71=1 

Here oi, . . . , ajv are distinct real eigenvalues, N < D, and Pi, ... , Pn are Her- 
mitean projection operators, with the properties that P m P„ = for m ^ n, 
and 

N 

i = E p « • ( 10 ) 

n=l 

This implies that 

N N 

(A) = Y,Pna n , var(A) = - ( A )) 2 - ( U ) 

n— 1 n— 1 

where p n — (P„) = Tr(pP„). According to the probability interpretation, the 
possible results of a measurement of A are the eigenvalues on, . . . , ojv, and p n is 
the probability of the result a n in the state p. 

The fact that the spectral theorem for complex Hermitean matrices is valid 
for all real symmetric matrices, with no more than the obvious changes in word- 
ing, allows us to generalize standard quantum mechanics by admitting as ob- 
servables all the real symmetric matrices. 



5 



3.2 States 
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In this generalization we have two options for choosing the set of states. The 
straightforward choice is to admit all real unit vectors as possible pure states of 
the system, and all real symmetric and positive definite matrices of unit trace 
as possible mixed states. This enlarges the set of possible states as compared to 
standard quantum mechanics, since not all such real density matrices are com- 
plex linear. In particular, it doubles the total number of states in the system, 
and it doubles the degeneracy of the spectrum of all standard observables, i.e. 
those that are complex linear. It seems that the degeneracy doubling is unphys- 
ical, at least if we want to describe systems that are well described by standard 
quantum theory. 

The second option, more interesting from the physical point of view, is to 
admit exactly the same states as in the complex theory. It means that we 
enlarge the class of observables, including all real symmetric matrices, but we 
admit only those density matrices that commute with J. For example, with J 
as in equation (0) the most general physical density matrix has the form 



(12) 



with 2(a +0) = 1, a > 0, j3 > 0, and a/3 - 7 2 - 8 2 > 0. 

This has the advantage that the physical degeneracies are unchanged, but 
it also has the somewhat strange consequence that there are no pure states in 
the theory. In fact, the pure states in the complex theory correspond to mixed 
states in the real theory. One way to see this is to observe that the meaning 
of the trace is different in the real theory as compared to the complex theory, 
because the number of basis vectors is doubled. Thus, if Tr p — 1 when the 
density matrix p is regarded as a complex D x D matrix, we have Tr p — 2 
when the same p is regarded as a real 2D x 2D matrix. Therefore the proper 
correspondence is that the complex density matrix p must correspond to the 
real density matrix p/2, which can never represent a pure state since it has no 
eigenvalues larger than 1/2. 

In this generalization of quantum mechanics as a theory defined on H 2D , 
every observable will have a complete set of 2D real eigenvalues and eigenvectors. 
However, one eigenvector alone does not represent a physical state. If we say 
that two physical states p and p' are orthogonal when pp' = p' p = 0, the 
maximum number of orthogonal physical states is D. Just by counting we see 
that if an observable A has more than D different eigenvalues, not every one of 
these can possess its own "physical eigenstate", in which a measurement of A 
gives this particular value with probability one. 

More explicitly stated, if an observable A does not commute with J, then it 
will have at least one eigenvalue which is not sharply realized in any physical 
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state, and conversely, there will exist no complete set of physical states such 
that var(A) = in every state belonging to the complete set. 



3.3 Poisson brackets 

To the classical Poisson bracket {A, B} of two classical observables A and B 
corresponds the "quantum Poisson bracket" 

{A,B) = ~[A,B\=~{AB-BA). (13) 

The proper way to write the same quantity in the real formulation is 

{A, B} = A£IB — BQ.A , (14) 

with ft = —J/h. The antisymmetry, {A, B} = —{B, A}, and the Jacobi identity, 

{A, {B, C}} + {B, {C, A}} + {C, {A, B}} = , (15) 

are easily verified. The most important property of the matrix £1 is that it 
is antisymmetric, Cl T = —SI, because that ensures that {A, B} is symmetric 
whenever A and B are both symmetric. Another way to write the relation 
{A, B} = C for symmetric matrices A, B and C is as 

[-JA,-JB]= -HJC . (16) 

This is then a commutation relation in the Lie algebra of the symplectic group. 

The real density matrix p must in general be explicitly time dependent and 
satisfy the Liouville equation 

£-| + t »H)-0. (.7) 

Here dp/dt is the absolute time derivative, dp/dt is the explicit time derivative, 
and H is the Hamiltonian. Thus the explicit time dependence of p is given by 
the equation of motion 

^ = {H, P } = Hnp- P nH . (is) 

The equation of motion must preserve Trp. A sufficient condition is that either 
H or p commute with fi, because then we have either pflH = pHfl or pflH = 
QpH, and in both cases 

= Tr(Hn P - pVLH) = . (19) 

Thus, if we accept all symmetric and positive definite matrices of unit trace as 
density matrices, we should impose the condition on the Hamiltonian H that it 
commute with J — —Ml. 
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We should impose the same condition on H even in the case where we accept 
only density matrices that commute with J. The point is that the equation of 
motion must preserve the condition of commutation with J, that is, the Poisson 
bracket {H, p] — HVlp — pflH must commute with J. A sufficient condition, 
when p commutes with J, is that H also commutes with J. 

When H commutes with J, and is not explicitly time dependent, the equa- 
tion of motion can be integrated explicitly to give 

p(t)=U(t)p(0)U(-t), (20) 

where U(t) — e~£ JH is the unitary time development operator. 

In conclusion, not every real symmetric matrix is an acceptable Hamiltonian 
in the generalized quantum mechanics as formulated here. It is necessary, or at 
least natural, to require the Hamiltonian to be a complex linear matrix. If we 
also require the density matrices to be complex linear matrices, it would seem 
that we are back to the point of departure, which was the standard complex 
quantum mechanics. However, we have enlarged the class of observables, even 
though we do not accept the new observables to be Hamiltonians governing the 
time development. 



4 The harmonic oscillator 

Time reversal was mentioned in the introduction as a motivation for the pro- 
posed generalization of the complex formalism. A more unconventional example 
of the generalizations that become possible, is the representation of the canon- 
ical Poisson bracket {x,p} — I in any finite and even dimension. For example, 
in the four dimensional case considered above, any two positive lengths £1,6 
define a representation of the form 
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The generalization to any even dimension 2D, or to infinite dimension, is ob- 
vious. Then both x and p anticommute with the imaginary unit J, as op- 
posed to standard quantum mechanics where they commute with J. With 
the most general physical density matrix p, equation (|l^), this gives (x) = 0, 
(Aa;) 2 = (x 2 ) = 2(a6 2 + /36. 2 )j an( l similarly for p. Thus the Heisenberg uncer- 
tainty relation holds in every physical state, 

AxAp = h^atf+K?) (|? + |?) 

= ^(« + /3) 2 +a/?(!-|) 2 >^. (22) 
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With these definitions the Hamiltonian of the harmonic oscillator of angular 
frequency u, 

v 2 1 

H = — + - muj 2 x 2 , (23) 
2m 2 v ' 

is diagonal and has the energy eigenvalues 

E < = *S? + 5 m " v ' < 24 > 

For any given value Ei > htu/2 this equation has the positive solutions 
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2Ei ± JAE 2 - h 2 iu 2 

N — ^ — ■ (25) 



By the obvious generalization, we may assign to the harmonic oscillator any 
finite number, or an infinite number, of arbitrary energy levels above the lower 
bound Tiu>/2. 

Note that the Hamiltonian H does commute with J, even though x and p 
here do not, implying that the time development operator e~^ JH is unitary. 
On the other hand, the operator e~^ Jp , representing translation in space by 
a distance d, is not unitary when d ^ 0, but only symplectic. The problem is 
that it is not orthogonal, because the matrix Jp in the exponent is symmetric 
rather than antisymmetric. Thus it does not conserve probabilities, and is not a 
symmetry transformation, as is evident from the fact that the position operator 
x has a discrete spectrum. 

The case £i = £2 = £ m equation (|2l]) is interesting because it is a realization 
of the canonical relation {x,p} — I, but is also, in a certain sense, a fermionic 
quantization of the harmonic oscillator, with 

h 2 

x 2 =t 2 , P 2 = J^> xp + px = 0. (26) 

In fact, the last relations are just one particular form of the canonical anticom- 
mutation relations 

aa T + a T a = I , a 2 = {a T ) 2 = . (27) 

To see this in more detail, let us introduce another antisymmetric matrix K so 
that it is an imaginary unit, K 2 = — J, and commutes with both x and p. Then 
wc take 

a=i*+itf P , ^ = ^x-iK P . (28) 



The Hamiltonian H as defined in equation ( |23| ) is just a constant, and a more 
interesting quantity is the usual Hamiltonian of the fcrmionic oscillator, which 
is 
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The time development operator with this alternative Hamiltonian is 
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In some respects this theory, where x and p anticommute with the imaginary 
unit J, is equivalent to another theory where the imaginary unit is K, which 
commutes with x and p. An explicit representation for K might be 
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with, for example, 
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If we define 



p = SpS- 1 = 
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then this is a density matrix which commutes with K instead of with J. If now 
p(t) = U'(t) p(0) U'(-t), then we have p(t) = U(t) p(0) U(-t), where 



u(t) = su'(t)s- 1 = 



(34) 



since SH' = H'S. The energy spectrum is the same in the two theories, for 
example we have 

Tr(pH') = Tt{pH') = 15hw . (35) 
However, the expectation values for x and p are not the same. For example, 

Ti(px) = , Tr{px) = 2(q - /?)£ . (36) 
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5 More than one degree of freedom 



With two degrees of freedom, referred to here by indices a and b, and belonging 
for example to two fields commuting with each other, the Hilbcrt space is a 
tensor product, 

H = H a ®H b . (37) 

However, the complex and real tensor products are mathematically different 
constructions. The complex tensor product of spaces of complex dimensions D a 
and Db has complex dimension D a Db and real dimension 2D a Db, whereas the 
real tensor product of spaces of real dimensions 2D a and 2Db has dimension 
4D a D b . 

The relation between the two types of tensor product can be understood as 
follows. In the complex case the following relations hold for tensor products of 
vectors, 

<j> <g> V = -i((i0) ® VO = -i(0 ® (WO) = -(WO ® (W) ■ (38) 

In the real case the four tensor products (f)®tp, ( Jcfi)®^, 4*®{J^) and {J(j))®{Jip) 
are linearly independent vectors, thus there exist two imaginary units J a = J® I 
and J b = I® J defined as operators on the real tensor product space Ti. It follows 
that H is a direct sum, H = H + ®H-, where H± = P±H 1 and the two operators 

P± = \{lTJ a Jb) (39) 

are complementary orthogonal projection operators. Each of the two subspaces 
TL± has dimension 2D a Db- Both J a and Jb commute with P±, and hence act 
within the subspaces Ti± separately. On H+ the relation J a — Jb holds, whereas 
J a = — Jb holds on H-. 

The complex tensor product space can be identified in a natural way with 
H+. Hence a physical density matrix p on the product space Ti must be a real 
(4D a Db) x (4D a Db) matrix that commutes with both J a and Jf,, and in addition 
it must have the property that p = P + p = pP + = P + pP + . 

If an operator x acts on H a and anticommutes with J, then the corresponding 
operator x a = x®I on H = Tia^Ti-b anticommutes with J a , but commutes with 
Jb- This means that x a maps H+ into 7i_, and vice versa. Such an operator 
could not be constructed at all within the complex quantum theory, because 
the subspace H- simply would not exist. It is essential for the construction 
that there exist two mutually commuting imaginary units J a and Jb on the 
real tensor product Hilbert space. The possibility of having operators that are 
antilinear with respect to J a and at the same time linear with respect to Jb, or 
vice versa, is at least a partial answer to one of the problems with antilinear 
field operators recognized by Stueckelberg and collaborators. 
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The operator x a above is "unphysical" in the sense that it maps the physi- 
cal space H+ out of itself. But similar unphysical operators are well known in 
physics. For example, a general isospin rotation does not respect the superse- 
lection rule for electric charge, it transforms a physical state into an unphysical 
superposition of states with different values of the charge. Another example is 
the relative position of two identical particles, which is an operator changing 
the symmetry properties of the two-particle wave functions. 

The generalization to more than two factors in the tensor product is easy, 
just introduce one new pair of projection operators for each new factor. If for 
example Ti = 7i a <8> T~tb ® W c , then define P± as above, and in addition 

Q± = \(lTJaJc) ■ (40) 

Since P e and Q v commute, for e = ± and r\ = ±, the products PJQ^ are also 
projection operators. Decompose H as a direct sum of subspaces Ji^ = P^Q^Ji, 
then in each such subspace we have J a = eJb = r/Jc- The physical subspace is 
it corresponds to the complex tensor product, since the relations J a = 
Jb = J c hold there. 
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